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Abstract 

Recently, it has been proposed that the Navier-Stokes equations 
and a relevant linear advection model have the same long-time sta- 
tistical properties, in particular, they have the same scaling expo- 
nents of their structure functions. This assertion has been investi- 
gate rigorously in the context of certain nonlinear deterministic phe- 
nomenological shell model, the Sabra shell model, of turbulence and its 
corresponding linear advection counterpart model. This relationship 
has been established through a "homotopy-like" coefficient A which 
bridges continuously between the two systems. That is, for A = 1 
one obtains the full nonlinear model, and the corresponding linear 
advection model is achieved for A = 0. In this paper, we investigate 
the validity of this assertion for certain stochastic phenomenological 
shell models of turbulence driven by an additive noise. We prove the 
continuous dependence of the solutions with respect to the parameter 
A. Moreover, we show the existence of a finite-dimensional random 
attractor for each value of A and establish the upper semicontinuity 
property of this random attractors, with respect to the parameter A. 
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This property is proved by a pathwise argument. Our study aims 
toward the development of basic results and techniques that may con- 
tribute to the understanding of the relation between the long-time 
statistical properties of the nonlinear and linear models. 
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linear advection models, random dynamical systems, stochastic analysis. 
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1 Introduction 
1.1 Motivation 

The GOY shell model [31] and [3H], and Sabra shell model [35] are some of 
the most interesting and most popular examples of simplified phenomenolog- 
ical models of turbulence. This is because, although departing from reality, 
they capture some essential statistical properties and features of turbulent 
flows, like the energy and the enstrophy cascade, and the power law decay of 
the structure functions in some range of wave numbers - the inertial range. 
We refer the reader to, e.g., [2], [Z], [H], [E], [2S|, and references therein for 
several descriptions and results. Often, in numerical or theoretical investi- 
gations, such models are driven by white noise forces. Both the stochastic 
GOY and Sabra shell models have the form 

dUn+ {l^klUn + bn{u,u)) dt = andPn, 71=1,2,... (l) 

where u„ (t) are complex valued, > is a parameter that represents the 
viscosity, fc„ = fco2" for some > are representing wave numbers, u{t) 
denotes the sequence (^))„>i, K (■, ■) is a complex valued bilinear func- 
tion of complex sequences u = {uj)-^^, that depends depending only on the 
variables Mn_2, tfn+i, w„+2 (where we impose the boundary conditions 
U-i {t) = uq (t) = 0). cr„ is a sequence of complex numbers, that are usu- 
ally chosen equal to zero for all n greater than some no (which describes the 
range of wavenumbers and consequently the length scales of external forces), 
{Pn)n>i is a sequence of independent complex valued Brownian motions. A 
rigorous theoretical analysis of the stochastic equation ([1]) and some of its 
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statistical properties have been investigated in [1], while other rigorous re- 
sults in the case of deterministic force have been developed in [H], [H], [U], 
[16]. The exact form of 6„ (■, ■) varies from one model to another. However, 
in all the various models in the sequel we assume that bn (■, ■) is chosen in 
such a way that 

oo 

(M,f)U„ = 0, (2) 

n=l 

for all square summable sequences u = ('Um)m>i ^ ~ ('^m)m>i • Equa- 
tion ([2]) implies a formal law of the conservation of energy in the inviscid 
(z/ = 0) and unforced form of ([1]). 

In analogy with the statistical theory of turbulence it is interesting to 
investigate the accompanying linear advection equation to equation ([1]), that 
is the linear auxiliary linear equation in the unknown w (t) = {wn (^))„>i 

dwn + {vklwn + hn [u, w)) dt = a„(i/3„, n = 1,2,... (3) 

where u is the solution of ([1]), and Wn (t) are complex valued functions. There 
is an extensive literature investigating the statistical properties of linear ad- 
vection (passive-scalar) equations in turbulent flows, which we do not pretend 
to cover in this contribution. We observe, however, that equation ([3]) is not 
the linearized version of equation of ([T]) about the solution u. This is because 
the term 6„ {w,u) is missing from ([3]), and an additive force still appears in 
the right-hand side of ([3]). Equation ([3]) should be considered as an auxiliary 
equation which, to some extent, may have similar statistical properties to 
those of equation ([1]), but is amenable to linear analysis (for instance the 
use of propagators). There is some numerical and heuristic evidence that 
some statistical properties of the solutions to equation ([3]), like the scaling 
exponents of the structure functions, are the same as those of the solutions to 
equation ([1]), see [1] and [5]. It is then of interest to understand the properties 
of the joint system 

dun + [vklun + hn {u, u)) dt = and(3n (4) 

dWn + [l^k^Wn + bn {u, w)) dt = andPn , 

for n = 1, 2, ... In addition, the following idea has been introduced first in [1] 
and proved rigorously later in [3]: one can symmetrize system (jl]) by means 
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of two additional terms as follows 

dUn + {l^klUn + hn (m, u) + A6„ {w , u)) dt = OndjSn (5) 
dWn + [vklwn + bn (m, w) + A6„ (w, lo)) = CTrtfi/?n , 

where A G M is a parameter, and to analyze the dependence on A of the 
properties of ([5]). For A = we recover (jlj). Observe that for A 7^ 0, setting 
V = Xw and multiplying the second equation by A, we have a perfectly sym- 
metric system for the pair {u,v), except for the force and initial conditions. 
Thus, to some extent, we would expect that u and Xw have similar statisti- 
cal properties for A 7^ 0. If we consider, for instance, the structure function 
Sp (kn) = {\un\^) (we do not specify at this heuristic level the meaning of the 
averaging procedure (.)), one might expect, as in the case of turbulent flows, 
that 

Sp {kn) ~ k~^^ 

for n lies in the so-called inertial range. The numbers Cp are called scaling 
exponents, which are universal in turbulent flows as the Reynolds number 
tends to infinity, i.e., as the viscosity tends to zero. Therefore, one possible 
definition of Cp is 

log5p(fc„) _ _ 

(;/,fc„)-^(0,oo) iOgkn 

where the limit is taken along a region of the form 1/°' < k~^ < u'^ for some 
a > (usually a = |). This is in order to ensure that the wavenumber 
considered are lying in the heart of the inertial range, as the viscosity tends 
to zero. Such statistical property, if it holds for w, it holds as well for Xw 
with the same value (p (and vice versa): indeed, if Sp^^ (kn) and Sp'^^^ (kn) 
are the structure functions of w and Xw, respectively, we have Sp'^^^ (kn) = 
XPSp""^ {kn) and lim(i,^fc„)^(o,oo) = 0, which imply the claim. Thus, if 
the scaling exponents Cp exist for both u and Xw (this assumption seems 
to be reasonable based on numerical finding in [1], [5] and the references 
therein) and are equal (which is reasonable to assume thanks to the symmetry 
u Xw described above), then they are equal for u and w. In summary, 
it is, therefore, reasonable to expect that some statistical properties like the 
existence and the value of scaling exponents, are the same for u and w, 
whenever A 7^ 0. 

Finally, it will be of great interest to show that such statistical properties 
depend continuously on A, as A ^ 0: if this is true, then the solutions of 
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([T]) and Q have the same statistical properties of the kind just described 
above. In particular, if this program is true, one is sure that results for the 
simpler linear model ([3]) can be translated to ([1]), which will be a remarkable 
breakthrough. 

1.2 Content of the paper 

The program above, outlined in [1] and [5], is composed of several steps, some 
of them are not easy to be justified rigorously. The first rigorous result has 
been obtained in [5] states that: in the case of deterministic forces, solutions 
of depend continuously on \ in C {[0,T]; H x H), for every given T > 0. 
Here H is the space of square summable sequences (t'n)„>i in C. This implies 
that the structure functions, defined as time average on any fixed finite time 
interval [0, T]: 



depend continuously on A. One of the limitations of this result of [5] that 
it considers deterministic forces. Here, we remove this restriction and prove 
the same result in the case described above of white noise forces. 

Several other issues have to be solved in order to be able to claim that 
the program described above is complete. One of the other major issues in 
[5] is that the statistics is being considered on finite intervals of time [0,T] 
instead of being considered on the attractor, i.e. as T — > oo. The existence 
of the limit as T oo, in the time average (definition of Sp{kn)) of the 
deterministically forced system is, therefore, one of these issues of [5]. We do 
not directly address this difficult problem here, in the stochastically forced 
case, but we content ourselves with a structural result about the infinite time 
horizon properties of ([5]): we prove existence of a finite-dimensional random 
attractor. This is a pathwise property, in the vein of the property of continu- 
ous dependence on A in C ([0, T] ; H) stated above for the deterministic case. 
We hope that this result, or the techniques involved in establishing it, may 
contribute to the understanding of the problem of the long-term behavior, 
i.e. T oo. Notice that we construct the random attractor for system 
for every A G M, hence if we take in particular A = the first component of 
the system is decoupled and thus the projection of the attractor on the first 
component is the random attractor of equation ([1]). Thus we prove in this 
paper the existence of a finite-dimensional random attractor for the stochas- 
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tic GOY and Sabra shell models, as a particular case of a more general result. 
However, the general result for system may help to prove further results 
on the relations between the statistics of the nonlinear and the linear cases. 

Due to the Ito nature of the previous equations, it is clear that other kind 
of analysis could be performed, in distribution and average sense instead of 
pathwise. This will be done elsewhere. We restrict ourselves here to purely 
pathwise properties. 

The paper is organized as follows. In section 2, we present the functional 
framework and prove pathwise well-posedness of system ([5]), and the continu- 
ous dependence of the solutions on A. In section 3, we give some preliminary 
results about random attractors and some of their properties. In section 4, 
we prove the existence of a random attractor for every coefficient A G M, its 
upper semicontinuity with respect to A; and finally that the random attractor 
has a finite Hausdorff dimension. 



2 Well-posedness and continuous dependence 
on A 

2.1 Functional setting 

Let us introduce the following spaces of complex valued sequences; we con- 
sider them as vector spaces on the field of real numbers. The space H is the 
space of P sequences over the field of complex numbers C: 

H = lu = {un)n>i '■ Un & C foT all 77, > 1 and \Unf < OO > . 

It is a Hilbert space with the inner product 

oo 

(u,v)jj := Re'^UnVr. 

n=l 



n=l 



and the norm given by = l^nl^- Let us recall that we have defined 

kn = 2"A;o, n > 1, with ko > given. We introduce now the Hilbert spaces 
D{A)cV CH defined as 

{oo 
Un\ < OO 
n=l 



6 



with norm \\u\\y = A;^ Moreover, for all a > 0, we define 

f oo 



oo 



n=l 



On the latter space we define the linear operator A" : D{A") G H ^ H as 

(A"ii)„ = C«n> for all ueD{A''). 
The operator A'^ is self-adjoint and strictly positive definite: 
{A'^u, u)h > kl"" \u\\ , for all u G 

We also observe that the inclusion maps of D{A) C V and V <Z H are 
compact embeddings. We finally introduce the bilinear operator B{-,-) : 
V X V —>■ H. For the GOY shell model it is defined as 

bn{u,v) := {B{u,v))^ 

A_ _ 1 /- - _ _ ^ 1_ _ 

:= ikn I -t^„-lU„+i - - {Un+lVn+2 + Vn+lUn+2> + 

For the Sabra shell model we define it as 

bn{u, V) := {B{U, V))^ : = = -fc^+l [(1 + S) Vn+lUn+2 + {2- 6) Un+lVn+2] 
+ -kn [(1 - 26)Un^iV„,+l - (1 + S)Vn-lUn+l] 
+ 3^n-l [(2 - S) Un-lVn-2 + (1 " 25) M„_2V„-l] 

(see [13], [H], [IS]), where 5 is a real number. In both shell models we impose 
the boundary conditions w_i = Uq = 0. What distinguishes the Sabra shell 
model from the GOY one is the dependence of the former on the parameter S, 
which is in charge for changing its character from the so-called 2d Turbulence 
regime to the 3d Turbulence regime, depending on the definiteness of the sign 
of a second (in addition to the energy) quadratic conserved quantity; see p] , 
[H, [El, [IE] and [35]. 

For both the GOY and Sabra shell models, the operator B (., .) is a bi- 
linear continuous operator from V x H to H, and also from H x V to H, a.s 
it will be stated in the next lemma. We also state its basic skew-symmetry 
property. 
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Lemma 1 There is a constant C > such that 

\B {u,v)\jj < C \\u\\y\v\jj , u&V,veH 

and 

\B {u,v)\jj < C \\v\\y \u\jj , v&V,u&H. 

Hence, B .) is a bilinear continuous operator from V x H to H, and from 
H xV to H. Moreover, 

{B{u,v),v)h^O 
for all u & V and v & H , or v & V and u & H. Equivalently, we have 

{B (u, v) , w)jj = - (-B (u, w) , v)jj 
for all u and v,w G H, or v,w & V and u & H . 
Proof. The first inequafity follows from the fact that 

oo /■ N oo 

^ kl \Un\^ \Vnf < ( sup \Unf ) ^ \Vnf < \\u\\y \v\% 
n=l \ n J 

and the second inequality follows similarly by interchanging u and v. Hav- 
ing proved these facts, the expressions in the last two identities arc all well 
defined. It is sufficient to prove the first identity, since it implies the second 
one because, ii u & V and v,w & H ox v,w & V and u G H, from the first 
identity and by the bilinearity of B we have 

— {B {u, V + w) , + w)jj 
= {B {u,v),w)jj + {B{u,w),v)jj , 

where we have used the fact that {B{u,v) ,v)jj — {B{u,w) ,w)jj — 0, by 
the first identity again. This implies the second one. One can also prove the 
converse. 

Finally, let us prove that 

oo 

J]Re[B {u,v)^vr^] = 0. 

n=l 
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For the GOY model we have 
— z 

— {B{u,v),v)jj 

Kq 

oo oo 

n=l n=l 

oo oo 

n=l n=l 

oo oo 

n=0 n=l 

oo oo 

- ^ 2""%„U„+iM„+2 + 2""%„M„+iTJ„+2 = 0. 
n=l n=— 1 

The computation for the Sabra model is very similar, see, e.g., [H]. The 
proof is complete. ■ 

We will also consider also the space V, the dual space of V, which can 
be identified as 

{oo 
u = (wn)„>i '■ Un & C for all n > 1, and A;^^ \unf < oo 
n=l 

with the norm \u\y, := Yl'^=i , u G V. It is clear that H C V, and 

V is the dual of V (with respect to H), with dual pairing between V and 
\^ defined as 



{u, v)y,y := Re ^ Wn'i'n, G V", w G V. 

n=l 



Observe that (m, v)jj = {u, v)y, y, when u E H, for every G 

It is easy to extend the operator ^4 as a bounded linear operator from V 
to v. One can also extend I? to a bilinear operator B {.,.): H x H ^ V . 
The definition is possible because 



oo / oo \ / 

\B{u,v)\y, = ^K^ \B{u,v)nf < C* ^ \Vn\'^ ( SUp \Un\ 
n=l \n=l J Vn>l 

<C. f kA ( kA = CMh\v\%- (6) 



in=l / \n=l 
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We also have 

{B {u,v),z)y,y = -{B {u,z),v)j^ 

for all u,v & H, z & V. Indeed, the identity is true for u & H, and v, z 
because in such a case {B {u,v) , z)y, y = {B {u,v) , z) jj and we may use 
Lemma [H Then we extend the result to G if by density of V in H. 
Define 

H = H X H, y = yxyand = x 

for a > 0. ^ ^ 

If X = (xi, xi) G H and y = {yi, 2/2) G H, we define the scalar product in 
H as 

<x,y >H=< Xl,yi >H + < X2,y2 >H 

and the norms in H and V as 

\^\% = I^iIh + \^2\h, X={xi,X2)eH 
ll^lll = ll^^lllv + 1^211^ ^ = {Xl,X2) G V. 

Moreover, define the linear operator A : D{A) G H ^ H, or also A : V ^ V, 
as Ax = {Axi, AX2) and, for every A G M, define the bilinear continuous 
operator Bx from V x H to H 01 from H x V to H a.s 

B\ (x.y) = {B {xi,yi) + XB (xs, 2/1) , S (xi, 1/2) + XB {x2,y2)) , 

where as usual we have used the notation x = (xi,X2), y = (2/1,2/2)- The 
main properties of the operator Bx are listed in the following lemma, whose 
proof is an easy consequence of Lemma [H 

Lemma 2 There is a constant C > such that 



Bx {u,v) 

and 

Bx {u, v) 



_< C 11^11^ , for every uEV,vEH. 



^ < C \\v\\rr \u\?r , for every vgV,ugH. 

H 



Moreover, 

Bx(u,v) ,v) ^ =0 



H 

for all u & V , V E H; or v E V,u E H ; also, 

Bx {u, v) ,w) ^ = - /Bx{u,w) ,v 

for all u eV and v,w E H, or v,w E V and u E H . 
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2.2 Well-posedness, stochastic flow and pathwise ver- 
sion in A 

Let (c7„) be a sequence of complex numbers such that 



for some £ > 0. This is a standing assumption for the sequel. 

Let Q be the space of continuous functions from M to H, null at zero, 
endowed with the metric of uniform convergence on compact sets. Let 
be the Borel a-field associated with Q. Denote by {W {t))^^^ the canoni- 
cal process defined on fl as W(t,u!) = uj{t), for every to & fl. Let P be 
a probability measure on (f2,jF) such that (VF(t))^>g and ^))t>o 
P— a.s. two independent Brownian motions in H with the same covariance. 
We call P a two sided Wiener probability measure and (W {t))^^^ a two sided 
Brownian motion. Such objects exist, for every given covariance operator, 
and play an important role in the theory of random dynamical systems, see 
[3j. Details on infinite-dimensional Brownian motions and their stochastic 
integration can be found in [20]. We will also denote by E the expectation 



For simplicity of the computations, and in analogy with equation ([T]), 
we assume that the components (VF„ {t))^^^, for all n > 1, of the two-sided 
Brownian motion have the form 



where (3n (t) are independent two-sided complex Brownian motions on {Q, JF, P) 
(with incremental covariance equal to one) and (cr„) is the sequence given 
above. 

On the probability space {Q, JF, P) consider the family of transformations 
{9t : I — > Q, t & M} defined as 9tUJ = uj {t + •)— a; (t), for every u & Q. They 
are measure preserving and ergodic with respect to P, and satisfy 9o = Id, 
9t+s = 9to 9s, for s, t G M, see [3]. 

Let (JF()jg]g be the filtration associated to (W^(O)teK i-^t generated by 
W{s) for all s<t). 

Given initial conditions uq,wo € H, let us first rewrite system ([5]) in the 
abstract form 




(7) 



on {n,j^, P). 



Wn (t) = anPn (t) 
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du^ = [-uAu^ - B{u^, u^) - XB{w^, u^)]dt + dW 
dw^ = [-vAw^ - B{u^, w^) - XB{w^, w^)]dt + dW 
u\0) = uo 

W^{0) = Wq 



We consider the above Cauchy problem on [0, oo). Using the notation of 
the previous section, system ([8]) can be rewritten as follows 

du^ + (^uAu^ + Bx {u^, u^) ^dt = dW (9) 

u\0) = Uo 

where (t) = {u^ (t) , (t)) , uq = {uo,wo), W (t) = {W{t),W{t)). As 
an introductory step, let us first give the usual definition of solution of ([9]); 
however, we will eventually need a more refined notion of solution, that we 
will introduce in Definition |4] below. 

Definition 3 Given uq G H, we say that a stochastic process u^{t,uj) is 
a solution of equation ^ if it is a continuous adapted process in H on 
JF, {J^t)t>Q 1 P) '^''T'd, for P-a. e. u E Q, 

u^{-,uj) G C{[0,T];H)nL'^{0,T;V) for all T > 

{u^{t, uj),i})g + V (u^{s, Lu),Ailj'j^ ds 

+ j (^Bx{u\s,uj),u\s,uj)) ,ij^_^ds 



{uo,^)H + (w{t,Uj),ij^^ 



for t > and ip ^ V . 



Notice that u^{s,uj) G V for a.e. s > 0, hence the integral of the bilinear 
term is well defined. 

The above definition is sufficient to analyze individual solutions, but the 
theory of random attractors requires the concept of stochastic flow: the P- 
negligible set where the properties of the above definition may not hold. 
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that is, it must be independent of Uq, and for P-a.e. a; G f2, moreover, we 
will also need continuity with respect to the intial value uq. In addition, 
here, we want to "vary" the parameter A independently of u: a priori this 
is not possible, again because the P-negligible set where the properties of 
the previous definition hold, may depend on A. Both problems can be solved 
because it is possible to perform a complete pathwise analysis of the equation. 
Let us, therefore, give a more appropriate definition of solution for which 
is relevant to the above mentioned issues. 

Definition 4 A stochastic flow depending on X eM, associated with equation 
is a family of mappings : H ^ H;t > 0,uj G A G m|, where 

E J-" is 6t-invariant and P {Q^) = 1, with the properties: 

1. for every A G M and Uq G H, {t,uj) t— > {p^{t,uj)uo (arbitrarily extended 
to all Lo E ri) is a continuous adapted process in H on (il, T , [Tt)^-^^ , P) 
and for every u E we have 

ip^{-,u)uo E C{[0,T];H)r]L\0,T;V), for allT > 0, 

and 



Jo ^ ' 



for t > and ip E V ; 

2. for every A G M and u E , <^'^{t, u) is a continuous map from H into 
itself, for all t > 0; and 

ip^{t + S,UJ) = (p^{t, OsUj) O ip^{s, LU) 

for all t, s > 0. 

To emphasize the role of in Definition |H we will consider stochastic 
flows depending on A G M, defined on the set Q^. 

We have the following result. The concept of uniqueness of stochastic 
flow depending on A means: if we have two stochastic flows, defined on two 
sets and then they coincide on a set E J-" such that P (fig) = 1. 
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Theorem 5 Under assumption there exists a unique stochastic flow de- 
pending on A G R, in the sense of Defintion^ associated with equation 



Proof. Step 1 (preliminary facts). Denote by e '^^^ the analytic semi- 
group generated by A (see, e.g., p6]). By the general theory of analytic 
semigroups or by explicit computation based on the spectral representation. 



for every a > we have 
over, notice that 



A^/^x 



< % for some constant C„ > 0. More- 



H 



H 



for all x & V. 



The process A'^/'^W (t) has if-components [A^/^W (t) , A^'/^W (t)) where 
A^/'^W (t) has complex components kf^anPn (t). Thanks to assumption ([7]) it 
follows that A'/^W (t) is an if-valued Brownian motion. Thus it is 7-Holder 
continuous, with respect to t, in H, for every exponent 7 < |, see [20]. This 
means that there exists a set ^l^^ G J-' such that P (fi^r) ~ -'- A^I'^W {t, co) 
is 7-Holder continuous for every exponent 7 < |, for every u G Q^y. The set 
is ^(-invariant, because Holder continuity is preserved by translation. 

Step 2 (auxiliary Stokes type problem). The pathwise analysis of equa- 
tion Qj requires a careful analysis of an auxiliary process. The process we 
are going to introduce is usually defined as 



but from this definition via a stochastic integral (which is a P-equivalence 
class) it is less easy to justify the ^t-invariance of certain properties, on a 
full measure set fl^. For this reason we adopt the following less intuitive 
definition. See [26] for further details on this approach. 

Let u G fl'^ be given throughout this step, where has been defined 
in step 1. The function t 1— > i'(t,a;) given by 



I{t, 



-vAU 



W{t,uj)+ uAe-'^'^'-'Uw{t,u;)-W{s,uj))ds, (10) 



is well defined and bounded in V, because (for e that is given in assumption 
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([7j)) we have 



e-^^'W (t, u) = ^i/2"^/2g-^At^./2 (w{t,iu)-W (0, iu) 



H 



< 

< c 



j^l/2-e/2 ^-uAt 
1 



H 



(w{t,Lu)-W{0,uj) 



H 



fl/2-e/2 



for every /3 < | and a suitable constant C > that depends on [3 and u. 
Observe that in the last estimate we used the details described in step 1 
above, in particular, the Holder continuity. Similarly 



j^^~vA{t-s) uj)-W (s, u)] ds 



V 



< 



< c- 



^3/2^e/2^-.Ait-s)^e/2 uj) - W {s, u) 

1 



ds 



H 



{t — s 



,3/2-e/2 



{t-sfds. 



The above estimates imply that is bounded on the interval [0,T], 

for all T > given, and the bound depends on T and u. 

With some additional minor effort one can show that the map t i— >■ 
uj) is continuous in V . We may write uj) componentwise: ^(t, uo) = 
, a;),z(^)(t,a;)) where, in the case when (T„ 7^ 0, 

zl:\t, u) = e'^'^'anPn {t, u) + f ukle""""'^'-'^ {a^Pn {t, uj) - anPn (s, uj)) ds 



with j3n{t,uj) defined as a^^Wn{t,uj) (if (t„ = 0, then Zn"* = 0). From the 
componentwise identity it is easy to deduce that 

{z{t,u),^)f,+ uh{s,u),Ai;)_ ds= (W{t,uj),ij)^ (11) 



vy 



H 



for allt > and i/j eV. 

Step 3 (auxiliary Navier-Stokes type random equation). Let u G 
be given, and let z{t,u)) satisfy (ITOl) or (ITT!) . Let us introduce the auxiliary 
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random differential equation 



^ ^J' + uAv^it, uj) + Bx{v^{t, u) + z{t, uj),v^{t, u) + z{t, u)) = (12) 

for t > 0. We say that v^{-,uj) is a weak solution of ffT^ if it belongs to 
C([0, r]; ^) n L2(0, T; V), for all T > 0, and if in addition it satisfies 

{v\t, uj),^)~+ V lv^{s,uj),M^^ ds 



+ 



for every t > and if) & V . 

For every uj G VL^^ and A G M, there exists a unique weak solution 
v^{-,uj) = v'^{-,u!,uo) of equation (fT2ll and it depends continuously, in 
C {[0 , T]; H) n L'^ {0 , T; V) norms, for any given T > 0, on the initial condition 
Uq in H. A full rigorous proof of this statement is very long, but at the same 
time it is very classical. Similar detailed proofs are given, for instance, in [1], 
[25], and in [12j or (3^ in the case of the classical Navier-Stokes equations 
(i.e., when z = 0). The rigorous detailed proof is based on the Galerkin 
approximation procedure and then passing to the limit using the appropri- 
ate compactness theorems. We omit these details which can be found in the 
above references. Instead, we present here the formal computations which 
lead to the basic a priori estimates, this is in order to stress the role played 
by 'z. Formally, if = v^{t,u) is a solution, then from various estimates 
and properties stated in Lemma [2] we have 

<C\v'\^\my{\v'\g+\z\H)- 

On a given interval [0,^], ||2; (■,a')||y and |z(-,co')|^ are bounded (see step 
2 above; they are bounded by a constant depending on a;), hence there is 
C (u) > such that 
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which imphes, by Gronwall lemma, a bound in terms of C (uj) and l^ollf 
for supjgjQ ~ and (t, a;) ||~ (it. These are the basic a priori 

bounds for the existence of weak solutions. For uniqueness and continu- 
ous dependence on initial data, we consider (■, and V2 {-ji^) two weak 
solutions, and we set {t, u) = Vi (t, u) — V2 (t, w), then formally we have 

+ uAy^ + {y\v^ + z)+Bx {v^ + z, y^) . 



dt 



, a; J I ^ on a 



Thus, by Lemma [2] and the boundedness of \z{-,Lo)\fj and 
given [0, T], we formally have 

<^l^1i?ll^iy(l^i|i? + l%) 

<2\\y\\v + -cMV\H 

for some constant C{v,uj) > 0. Again by Gronwall Lemma, uniqueness and 
continuous dependence on initial value follow. 

Step 4 (existence of the stochastic flow). For every 
and Uq & H, define 

Uj)uo = {t, UJ, Uq) + z{t, Uj) 

where ■u'^ (■, w, Mq) is the unique weak solution given in step 3 and z{-,ljj) is 
defined in step 2. The set Q^r is ^t-invariant and P = 1. Property 1 of 
Definition mis a direct consequence of the analogous properties ofv^ (■, f^, Mq) 
and z{-,uj) proved in steps 2 and 3. As to property 2 of Definition UJ given 
A G M, c<j G Q^, t > 0, the continuity of (/?'^(t,t<;) in H is a consequence of 
the continuous dependence of (■,a;,'Uo) on uq, see step 3. The property 

ip^{t + s, uj)uo = ip^{t, 6sU!)lp^{s, uj)uo (13) 

for all t,s > follows from the uniqueness statement of step 3. In order to 
prove this claim, let us write, for a given s > 0, the equation satisfied by 
the two functions t h-» ifi^(t + s, uj)uo and t h-* if^{t, 9sUj)ip^{s, uj)uq for t > 0. 
We know that Lp^{t,u)uo satisfies the weak form of the equation given in 
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Definition SI From it, for tlie function y (t) := ip^{t + s,uj)uo, we have 



t+s 



/n \ / V.V 



dr 



+ 



t+s 



Bx {v^{r, uj)uo, ip^{r, uj)uo) ,ip\dr 

I M 



{uo,ij)^ + {W{t + s,u),ijJ^ 



for all t > and ifj eV. Hence, 



i^{yir),A^p)^ dr+ (Bxiyir),y{r)),^lj)dr 

Jo ^ ' ^ y Jo ^ ' ^ 

That is, 

{yit),i')H+l (y (r) , Atp) _ dr + I lBx{y{r),y{r)),%l))^dr 



H 



vy 



H 



H 



Recall that W {t, u) = {W {t, u) , W {t, u)), W {t, u) = u (t), and 

W {t, OsUj) = OsUJ (t) = uj{s + t) -u (s) 
= W{s + t,u)-W{s,uj). 

Therefore, y satisfies 

{yi.'t) ,'^)h+ I (y (r) , Alp) dr + [Bxiy (r) ,y (r)) dr 



vy 



H 



H 



This is the same equation satisfied by the map t ^ ip'^{t,6sUj)(f{s,uj)uQ for 
t > 0. Since this equation corresponds to the auxiliary equation of step 3 
through the transformation via z{t,uj) (the detailed argument is the same 
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as the one given below in step 5, but easier since here = ^l^/, see below) 
and we have uniqueness for the latter, thus we also have uniqueness for the 
former. This proves f|T3|) . 

Step 5 (uniqueness of the stochastic flow). Let Lp^{t,u!) be a stochastic 
flow depending on A G M, associated with equation ([9]), defined on a 9t- 
invariant full measure set fi^- Consider the ^(-invariant full measure set ^^/^ 
described in step 2 above. The set Q^HQ^^ is 6'rinvariant and P {Qi fl ~ 
1. For all a; G n define 

{t, LO, Mo) := ^i{t, uj)uo - z{t, uj) . 

From the properties of tp'^ and z it is trivial to check that (-jWjMo) is a 
weak solution of the auxiliary equation of step 3. Of course we have 

Now, let ip2{t,u!) be another stochastic flow depending on A G M, as- 
sociated with equation ([9]), with its 6'f-invariant full measure set ■ The 
function 

V2 (t, UJ, Mo) '■= V^2(^) ^)Uo — 'z{t, Uj) 

defined for a; G ^ weak solution of the auxiliary equation of step 

3. Thus, for u E H n9, H Q'^ we have 

because of the uniqueness of solutions for equation f|T2|) . for every tu G fl^r. 
Thus 

for all G fi? n f2° n Since P n n fi^) = 1, the proof is complete. 
■ 

In step 5 of the previous proof we have obtained also the following repre- 
sentation result, which will be useful in the next section. 

Corollary 6 Let ip^{t, u) be a stochastic flow depending on X eM., associated 
with equation defined on a Ofinvariant full measure set On the 9t- 
invariant full measure set Q^r described in step 2 of Theorem one can 
define the functions z{t,uj) and {t,u,UQ) according to steps 2 and 3 of 
that proof. Then, on the Ofinvariant full measure set := fl we 
have 

ip^{t, Uj)uo = {t, UJ, Uq) + z{t, u) . 
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2.3 Continuous dependence with respect to the pa- 
rameter A 

As above, we assume condition ([7]). The uniformity in the initial condition 
of the next statement will be used to prove the upper semicontinuity of the 
random attractor with respect to the parameter A. 

Theorem 7 Let (fi^{t,uj) he the stochastic flow that was established in The- 
orem and Corollary 0, associated with equation ^ and depending on the 
parameter A G M. Let E J^, P = 1, be a 6t-invariant set where all the 
properties of Definition^ and CorollarylB hold true. Then, for every uj e 0°; 
we have 

lim sup sup \Lp'^{t,uj)uQ — ip'^°{t,uj)uQ\^ = Q 

A — ►Ao S„gBO< t<T 

for all T > 0, Ao G M and all bounded sets B d H . 

Proof. We prove the theorem only in the case Aq = 0, the general 
case being the same. The elements w G and T > are given and fixed 
throughout the proof, as well as the bounded set B C H. 

Step 1 (preparation). Denote by [u^{t,Ljj,uo),w^{t,uj,uo)) the decom- 
position of (p^{t,uj)uQ in H = H X H, and by {uq,Wo) the decomposition of 
the initial value uq- Where it is necessary, we will shorten the notation and 
write ('u'^(t), w^{t)^ and apply analogous change of notation to other similar 
quantities. 

From the weak integral equation in Definition H] and the definitions of A 
and Bx we have 

Jo 

+ [ {B {u\s),u\s)) + XB {w\s),u\s)) ,^,)^ds 
Jo 

= {UO, i>l)H + (^(^' ^^)h ' 

for all t > and ifji eV and 

{w^{t),ij2)H+ I l^{w^{s),Aij2)yy,ds 

Jo 

+ / {B{u^{s),w^{s)) +XB{w^{s),w^{s)) ,^2)Hds 
Jo 
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for all t > and '02 ^ Let us define the new function 

q^{t, LU, uq) := u^{t, u, Uq) + Xw^{t, Uo) 
and the corresponding difference 

p^{t,uj,Uo) := q^{t,uj,Uo) - q'^{t,uj,Uo) = q^{t,LO,Uo) - u'^{t,uj,Uo). 
The above quantities are solutions, respectively, of 

{q\t),^i)^+ [ u{q\s),A^,)^^,ds 
Jo 

+ [ {B{q\s),q'^is)),ij,)^ds 
Jo 

= {uo + Xwo, A)h + (1 + ^) (^(^' ^).^i)h > 
for t > and ipi G V , and 

(p^(t),V'2>^+ / u{p\s),Aij2)yy,dS 

Jo 

+ / {B{q\s),p\s))+B{p\s),q\s))-B{p\s),p\s)),^lj2)^ds 
Jo 

for t > and 1^2 e V. 

Step 2 (bound on q^). Let us prove next that 

sup sup sup \q^{t, u, uo)\rr < 00. 
Ae[-i,i] uoeB o<t<r 

Define 

V^{t, UJ, Uq) := q^{t, iO, Uq) - (1 + A) z{t, u) , 

where z{t, u) is any one of the two equal components of 'z{t, u) given in 
Corollary [6l We have 

{v\t),^i)^+ ! v{v^{s),A4j^)^^,ds 
Jo 

+ / {B{v^{s) + {l + \)z{s),v^{s) + {l + \)z{s)) ,i)i)^ds 
Jo 

= {uo + Xwo,'iJji)h , 
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for all t > and ifji G V. Formally, this implies 

Ijt l^'l^ + ^ ll^'ll'v- - I <^ + + ^) + + ^) ' 
and thus 

1 I \ , i2 /"* II \ . . ii2 . 1 



2 1^"^ W|h I i^)\\vds < ^\uo + Xwo\ 



H 



+ / \{B{v\s) + {l + \)z{s),v''{s) + {l + \)z{s)),v\s))^\ds. 
Jo 

Rigorously, the above inequality can be proved either by general abstract 
theorems (see [37]) or by taking finite-dimensional (i.e. with finite many 
components) test functions ipi, performing the computations at the finite 
dimensional level and then taking the limit, which can be justified because 
the map 



s ^ {B {v^ (s) + (1 + A) ^ (s) y (s) + il + X)z (s)) y (s))^ 

is integrable. Thus, from Lemma [2] and the bounds on z {t,uj), given in step 
2 of the proof of Theorem [5l we have, for t G [0,T], 

<2 / \{B{v^{s) + {l + X)z{s),z{s))y{s))^\ds 
Jo 

<C f \v' {s)\^ is) + il + X)z {s)\^ \\zis)\\yds 
Jo 

Ciuj) [ \v\s)\^{\v'{s)\^ + C{u;))ds, 
Jo 



< 



where C (cu) depends on the bounds of the relevant norms of 2 (-, u;) over the 
interval [0, T], which in principle is also depending on T. By Gronwall lemma 
we deduce 

{t,u,Uo)\^ <C{uj)- (0,a;,no)|^ 
on [0, T], for a new constant C (uj). This implies 



\q^ {t,uj,Uo)\^ < 2 \z{t,u)\^ + C{uj)- \uo + Xwo\ 
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H ' 



on [0,T], and the claim of this step is proved, using again the bounds on 
z{-,uj) over the interval [0, T]. 

Step 3 (convergence of p^). Next we prove that 

lim sup sup \p^{t,uj,Uo)\ =0. 

In one sentence, this is a consequence of the various bounds, that we have 
established previously, and the fact that the initial condition Xwo and the 
forcing term XW{t,u) converge to zero, as A — > 0. Define 

v^{t,uj,Uo) := p^{t,uJ,Uo) — Xz{t,uj), 

it satisfies 

Jo 

+ [ {B {q\s),p\s)) + B {p\s),q\s)) - B {p\s),p\s)) ,^2)^ds 
Jo 

for alH > and 7/^2 G V . By virtue of Lemma [1] and Lemma El we have 

\{B {q\p') ,v')^\ = \{B {q\ Xz) ,v')^\ < C\X\ \\z\\y \q'\^ \v'\^ . 

Thus, using both the bound on 2; (-, c<j) on [0, T] and the bound of the previous 
step, there is constant Ci (u) such that 

\{B{q\p'),v')J<X'C^{co) + \v'\l, 
for all Uq E B and t G [0, T]. Similarly, 

1(5 (p^g^),^;%| = 1(5 [v^ + Xz,v^),q^)^\ 

<C\\v'\\y\q\{\v\ + \X\\z\^) 

for some constant C2 (oj), and 

\{B {p\p^) ,v^)^\ = \{B {v' + Xz, Xz) ,v')J 

<C\X\ ll^ll^l^^l^d^^l^ + IAII^IJ 
<|A|C3M|^^|^(|^^|^ + |A|C3H) 
<\X\Cs{u;)\v't + \X\'Cs{uj) , 
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for some constant C3 (cj). Hence, from the equation in weak form for (and 
similarly to the proof of step 2 above) we deduce, for A e [—1, 1], 



\v^{s)\\ydS 



for some constant C4 (oj). By Gronwall lemma we get 

lim sup sup a;, Uo) I IT = , 

^ — *OSoeBO<t<T 

which implies the claim of this step. 

Step 4: (convergence of w'^) . With the notation (^'*'(t, a;, So) := w'*'(t, a;, wq)- 
u'°(t, a;, Mo), let us prove that 

lim sup sup |^^(t, a;, uq) I „ = 0. 

^ — '^uoeBa<t<T 

We have 



[B {u\s),w'is)),i;)^ds 

- f {B {u\s),w\s)) + XB {w\s),w\s)) ,'ilj)^ds 
Jo 

= [ {B{q%s),w%s)),ij)^ds- f {B{q\s),w\s)),i;)^ds 
Jo Jo 

= -/ {B{p\s),w%s)),^)^ds- f {B{q\s),e{s)),i^)„ds, 
Jo Jo 
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for all t > and -0 G V^. As in the previous steps we deduce that 

l\C"it)fn + '^j\\^"is)\fvds 

< f\{B{p\s),w\s)),e{s))^\ds 
Jo 

Jo 

<C f\w'{s)\^\\e{s)\\y\p\s)\^ds. 
Jo 

Hence, 

This implies the claim of this step. 

Step 5: (convergence of u^). We simply notice that 

\u\t)-u%t)\^=\q\t)-q'it)-\w\t)\^ , 

therefore, by the results of steps 3 and 4, we have 

lim sup sup \u^{t,u,uo) — u'^{t,uj,uo) \ „ = 0. 
^ — *OSoe-BO<i<r 

The proof is complete. ■ 

3 Random dynamical systems 

In this section we recall few definitions from the theory of random dynamical 
systems. For general notions and results see [3], and see for analogous 
concept for non-autonomous dynamical systems. Here we mainly refer to 
specific notions from [19]. 

3.1 The basic set-up 

Recall from the previous section, the following notation: let {Q, J-", P) be a 
probability space and {6t : i — > Q, t ^ R} a family of measure preserving 
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transformations such that {t,uj) i — > Ot{uj) is measurable, = Id, 6t+s = 
6t o 9s, for s,t G M. The flow {9t)tm together with the probability space 
{Vt, T , P) is called a measurable dynamical system. Furthermore, we suppose 
that the map 9t is ergodic. 

Definition 8 Let {X,d) be a Polish space (i.e. complete separable metric 
space) and B its Borel a-algebra. Let = [0, oo). A map 

ip:R+ xnx X I — > X 

{t, U, x) I > ip{t, U!)x 

is called a measurable random dynamical system (RDS) on X over [Q, T , P, 9t) 
if the following properties are satisfied 

1. if is (i3(]R+) -B, -B) measurable, where i3(R+) is the Borel a- 

algebra on R+; 

2. 

(p{t + s,u) = (p{t,9siu) o ip(s,uj), (14) 

for all t, s G M"*" and Lp{0,uj) = Id, for all u & Q. Property [l4\ ) is 
called the Cocycle property. 

An RDS (f is said to be continuous or differentiable if for every fixed 
{t, oj) G xfi, Lp[t, uj) : X I — y X is continuous or differentiable respectively. 

Instead of assuming f|T^ for all u; G fi, it suffices to assume it for all u 
from a measurable ^(-invariant subset of full measure. 

3.2 Attraction, absorption and invariance 

Let {X, d) be a metric space. For two nonempty sets A, B d X , v^e recall 
the Hausdorf semi- metric du {A, B) = sup^g^ infygs d (x, y). 

We observe that dn restricted to the family of all nonempty closed subsets 
of X is a metric, see [9]. 

Definition 9 Let be a measurable space and let {X,d) be a Polish 

space. A set valued map K : VL — > 2^ , taking values in the closed subsets 
of X, is said to be measurable if for each fixed x E X , the map uj i — > 
dnix, K{u)) is measurable. The map K is often called a closed random 
set. 
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Definition 10 Let ip : M"*" x f2 x X such that {t,uj,x) i — > ip(t,u)x G X 
be a measurable RDS on a Polish space {X, d) over a measurable dynamical 
system P, Ot). A closed random set K is called (/^-forward invariant 

if for all G fi, 

'{){t,uj)K{uj) C K^Otto) for all t > . (15) 
A closed random set K is called strictly (/^-forward invariant if for all 

ip{t,uj)K{uj) = K{9tuj) for all t > 0. (16) 

Remark 11 By substituting 9_tUJ foru in Definition lT^ we get the following 
equivalent version of Definition ITU A closed random set K is called strictly 
99-forward invariant if for all for all uj & Q, 

ip{t, e.tio)K{e^tio) C K{uj) for all t > (17) 

or respectively 

ip{t, e^tUj)K{e_tUj) = K{uj) for all t > 0. (18) 

Definition 12 A closed random set K is said to absorb the set B d X , B 
is fixed non-random, if there exists a random variable tsi^j) such that, for 
P-a.e. oj dVt, 

ip{t, e^tUj)B C K{uj) for all t > tB{uj). (19) 

The smallest tsiuj) > for which / fTPj) holds is called the random absorption 
time of B by K. 

Remark 13 Note that ip{t, 6-t^)x can be thought of as the position of the 
trajectory at time 0, which was in x at time —t. 

Definition 14 For a given closed random set K , the uj-\\m\t set of K is 

defined to be the random set 

^k{uj) = n U <^-t^)K{9^tUj). (20) 

n>0 t>n 

Remark 15 1. A priori Kxi'^) can be an empty set. 
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2. We have the following equivalent version of Definition 



Ak{uj) = \ y e X : 3tn^ oo, {xn} C K{e^t„uj), lim ip{tn,0_t„uj)xn = 

(. n—^oo 



3. Since [J^~^^(p{t,9^t^^)K{9^t'^) is closed, then A/^(a;) is closed as well. 

Definition 16 A random set A{uj) is called a random attractor associated 
with the random dynamical system ip if, for P-a.e. uj & Vt, the following is 
satisfied: 

1. A{(jj) is a nonempty compact subset of X , 

2. (p{t,uj)A{Lo) = AiOtiu), y t>0, 

3. for every B G X bounded (and non-random) 

lim dHivit,9^t^^)B,A{uj)) = 

t >oo 

The following theorem about the existence of random attractors is due 
to Crauel and Flandoli [T^ . 

Theorem 17 Suppose there exists a closed random set D which is absorbing 
every bounded non-random set B G X , and for which D[uj) is a compact 
subset of X for P-a.e. u gVL. Then, the set 

A{uj) = U Kb{oj) , 

BCX 

is a random attractor for ip. Where the union above is taken over all the 
bounded and non-random B G X, and Ab{uj) is the u-limit set of B. 

Remark 18 In Crauel [17] it is shown that, under the ergodicity assumption 
on 6t, there exists a compact set K{uj) G X such that, for P-a.e. u G Q, the 
random attractor is the uj— limit set of K{uj), that is, 

n>0 t>n 
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3.3 Random attractor dimensionality 

We are interested in the property of finite Hausdorff dimensionality of the 
random attractor. The two most relevant results for this purpose are the 
works of Debussche p2] , jJSj . We apply the result of the second one of these 
papers, based on a property called random squeezing property, which was 
inspired by the squeezing property in the deterministic case that was intro- 
duced in [28] (see, also, [13] and [M])- This property is also used in the proof 
of finite number of determining modes. The fact that the random attractor is 
not uniformly bounded makes the corresponding random squeezing property 
depend exponentially on a random variable. However, an ergodic argument 
will make it possible to work with this weaker property. 

Definition 19 (fEEl) Let H be a separable Hilbert space with norm |.|^; 
(p{t,u) a random dynamical system in H with random attractor A{u). We 
say that (f{t, u) satisfies the random squeezing property if there exist a random 
variable C^i^uj), a finite- dimensional projector U in H, and positive numbers 
/i, 6 such that, for P-a.e. u & Q, 

\Uip{t,Lu)uo-Ilip{t,uj)vo\H<e^o0.iOs^)ds\^^_^^\^ (21) 

and 

\{I -U){if{t,u;)uo-^{t,u;)vo)\H < (e"^* + c^e^o - t;o|/f (22) 

for every t > 0, and every Uq, Vq G A{uj). 

Theorem 20 (122^ ) There exist absolute constants Ki, K2, such that if 
(pitjUj) is a random dynamical system that satisfies: 

(i) the random squeezing property, mentioned in Definition [7P|. with a 
random variable C^{u), a finite- dimensional projector 11 and two positive 
numbers /i, 5, 

and 

(a) the expected value with respect to the measure P 

E{C^{uj)) < 00, 5<Ki andfi> K2E {C^{uj)) , 

then, for P-a.s. uo & Vt, the random attractor A{uj) of (f{t,u) has finite 
Hausdorff dimension which is less than KsR{Il) log/2(n), where -R(n) is the 
rank of the projector U. 
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4 Application to the Shell model 



In Theorem |5] we have constructed, for every A G M, a random dynamical 
system (p^(t,u) associated with equation ([9]). In this section we prove the 
existence of the random attractor associated with (p''^(t, uj). At the end of the 
section we also prove that the random attractor, as a function of A, is upper 
semi-continuous. For the upper semi-continuity of deterministic attractors 
with respect to a parameter see, e.g., [32j . 

4.1 Auxiliary problem 

As in step 2 of section 2.2, we will introduce an auxiliary Ornstein-Uhlenbeck 
process. This is a slightly different process but will have mainly the same 
properties as the process introduced in step 2 of section 2.2. Following the 
steps of section 2, let G introduced in section 2. Let a > be an 
arbitrary constant. For to G let 



z''{t,tu) = / {uA + aI)e-^''^+"^^^'-'\W{t,u) -Wis,u))ds (23) 



which is well defined and bounded in V. 

The process 2^"(t), t G M is a Gaussian, stationary and ergodic process. 
Moreover, it is a solution of the equation 




— oo 



dl'^{t) = {-uA - a)z''{t)dt + dWit) 



(24) 



i.e. for all t G M and P— a.s. 




(25) 



In particular, for each component one has 




(26) 



Moreover, 
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n=l 

n=l V^-oo / 

°° / ft _ \ ■ 

«=1 \J-co / 



n=l 
oo „t 



1 2(z/A;2 + a 

n=l ^ "■ 



Furthermore, tends to 0, when a —>■ oo. In particular, there 

exists a^{i') > such that 



^ll^"Wllv'<^ for all a>a,, (28) 
where C^, is the constant in the inequality (Q. 

4.2 Absorbing compact set 

Let u G fl^r be given and let us introduce the random differential equation 



^^ + .Avit,u)= (29) 
at 

- Bxiv{t, u) + ^"(t, uj),v{t, u) + ^"(t, a;)) + ^^"(t, a;). 

Notice that in fact v depends on a, because depends on a. It is not difficult 
to prove, using a Galerkin method, that for each u G and to £ with 
VtQ^u) E H is given, there exists a unique solution v{t,uj) defined on [to, oo) 
to dSni) such that 

v{to,uj) =Vtoiuj) , (30) 

and such that 
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v{-,u) e C([to, oo), fl L\{t^, oo), y). 

We refer to [22] for more detailed computations. 
Let us define 

V?^(t, Uj)Uto = V (t, u) + 2" {t, LO) 
where v is the solution of (129|) with 

We would like to prove the existence of a compact absorbing set in H at 
time t = 0. Through this section we will take i? to be a bounded set in H, 
and that for any to £ we will assume u{to) G B; moreover, v is the solution 
of and dnni) with 

vtoiuj) = u{to,uj) - ^"{to,^). 
Let to < -1 and t G [-1,0]. 

Lemma 21 Let Vt^ G H and v he a solution of (1291) associated with the 
initial condition VtQ- Then, for all to < —1 and for all t G [—1, 0] 



J to 

(31) 

where 

fit) ■.= -^\nt)\%+^mt)\%, 

and C* is the constant in the inequality 

Proof. As before, the proof is formal and can be made rigorous by ap- 
plying the Galerkin approximation procedure. Let us take the inner produce 
in H of equation fl29l) with v to obtain 



(J/V ' — ' — 

<—,?;>+//< Axi,v>= - < Bx(v + z°,v + z°),v > +a <z",v> 
dt 
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Using Lemma [21 inequality and Young's inequality, we estimate the right- 
hand side of the above and get 



< Bx{v + z",v + z",v > \ = \ < Bx{v + z°',v),v > + < Bx{v + z°', z",v > \ 

= + > I = I - Bxiv + z'',v),z'' > 

= \-Bx{v,v),z'' > -Bx{T,v),z'^ > I 

<C4\v\l\\r\\y+\\v\\y\I%) 



<CMW''\\v + ^\\vry + -jf\^%- (32) 

For the other term we apply the Cauchy-Schwarz and Young's inequalities 
we have 

|a<F,^T> I < alz^lfrlvlf. < ^\v\l + ^\z%. (33) 
Hence, we get that 

2 * + 2 IMII < C.Klk-ll. + + -fHl + (34) 

Using the Poincare-like inequality z^HvHl, > z^/cqI^I/^, we get 

We integrate over {tQ,t) and get 



I^^WI1^< l?^(to)ll, + ^* (2a||."(r)||^-^) \v{r)\ldr + j\r)dr. 

Using Gronwall Lemma we get the result. ■ 

Lemma 22 Let a > a*(i^), such that /[^) holds. Suppose that the assump- 
tions of Lemma IM\ are satisfied, and let f be given by Lemma IMl Let 

R{iuj):=l+ f /(.)e^^'(2^*i'^"('-)ilv-^)'^^ds 

J ~oo 
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then, the ball B{0, Ri{uj)) G H is an absorbing set at time t G [—1,0], for 
the system /[^) . 

Notice that in principle R\{uj) depends on a, however, if we fix a = ao : = 
2a^,(z/) (see (!28l) ) then R\{uj) will depend only on u and the statement will 
still be valid. 

Proof. Using the ergodicity of the process we have that 

z''{s)\\yds ^i„^_oo E\\z''{0)\\y, P - a.s. , 
Hence, there exists Sq{uj) < such that for every t^ < Sq{ijj), 

-r fw{s)\\yds<2E\\Tmv- (35) 

*0 J to 

Then, 




exp [t, + ^ j\cA\T{s)\\yd.^^ < exp [t, - 2ai?||?"(0)||^j j . 

Moreover, thanks to (I27|) and we have for all a > a* 

E\\z-mv<^^- (36) 
Hence, for every to < Sq{<^) and a > a* one has 



exp {t, + 1. j\c4z"{s)\\yds^^ < exp (^t^- 



In addition, see [25] for more details, one can easily prove that there exists 
an a.s. finite random constant Ct{uj) such that 

\z'"{t)\H < C7{u)\t\, for all t < -1. (37) 

Let us assume that the bounded ball B G H is inside a ball of radius pi, 
then 

\u(to)\fi — Pi) ^0 < 0. 

Hence, for every < So{uj) and a > a* 
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I- — / \ \ '2 1-^/ / \ I - — rv / 1 2 / koJy 

\v{t)\rr < \u{to)\fyexp [to-— ] + \z {tolfyexp [to 



J to 

+ 



8C* ^ 

Now choose Si{u) < such that (pj + Cj{u)\to\'^) exp (^^o^) < 1, for all 
to < si{u), and let us denote by tB{uj) = mm {so{uj) , si{uj)} , then we get 
that the integral inside R\{uj) is a.s. convergent and that for every t < tsiuj) 

mil < R\iu^)- 

Hence, the ball B{0, R\{u)) is an absorbing ball at time t. This completes 
the proof. ■ 

Lemma 23 Suppose that the assumptions of Lemma l2J\ are satisfied, and 
assume that a > and |so(a;)| is large enough such that (EQi o,nd l[35\) hold, 
respectively. In addition, assume that Vt^ G V , then for every t^ < tsiuj), 
there exists R2{uj), P- a.s. bounded, such that 



\\v{t)\\ldt<R2{iu), 



1 



where 



R^iu) := ^Rf\uj) + f_J{t)dt + j\{{u) [2C4z'^{s)\\y^^^ ds. 

Notice again that in principle R2{uj) depends on a, however, if we fix 
a = ao := 2a^,(z/) then R2{oj) will depend only on v and the statement of 
Lemma [23] will still be valid. 

Proof. Integrate fl34p over (—1,0), then use Lemma [2T] to estimate 
\v{-l)\\. ■ 
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Lemma 24 Assume the assumptions of Lemma are satisfied, then there 
exists R^luj), P-a.s. finite such that 

||^j(o)|||</?3H, 

where 



Rsico) := (^R,{u) + l'^(^2C'ARl{u;) + \I-{s)\l)\\z^{s 
X exp(^° 2C^XRI{lo) + . 

Proof. Let us take the inner product in H of the equation fl29l) with Aw, 
we get 



1 d 



Ivfy + iy\Av\g = - < Bx{v + z°',v + z°'),Av > +a < z°',Av > 



^' '~|2 I r~i2\~ I ~a\2 ||~ I ~Q||2 I 4a |~o|2 



H ■ 



In the above estimate, we have used Lemma [2] and the Young's inequahty. 
Hence, 

ft 

|2 ^ ll~/'„M|2 I / 0/^2|~/^\ I ~a/\\2 ||~/„M|2 



J s 

+ (^Cl\v{r)+T{r)\l\\T{r)\\l + ^|F(r)|| ) dr. 
Using Gronwall lemma, we get 

+ f (2Cl\v{r) + z-{r)\l\\T{r)\\l + —^^ 



V 

Therefore, for t = 0, we have 
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Il~^n\l|2 ^ M|2 r°2Cj|iJ(r)+2"(r)|2 dr 

I /"° /or^2|~/ \ I ~OL( m2 ||~a/ M|2 , '^^^ \~l m2 ^ 2CJ|^^(r')+^°(r')P- dr' , 

Now, we integrate over (—1,0) to obtain 



IRO)lH < (£11^^(^)111^5 + £(2C,^|?^(.) + z(.)1|||z"(.)||| + ^|F(.)||)rf.) 

xexp {J\cl\v{s)^z'^{s)\\d^ . 

Consequently, we use the estimate of the preceding lemma to complete the 
proof. ■ 

Lemma 25 Let (fi^{t,u) be a stochastic flow associated to equation de- 
fined on a Ofinvariant full measure set On the Ot-invariant full measure 
set Vt^r described previously, one can define 'z'^it.uj) and v{t,u,Uo), say for 
a given fixed a = := 2a^(z/) (see (d^jj- On the Ofinvariant full measure 
set n° = fi? n we have 

ip^{t,u)uo = v{t,u,Uo) +z°'{t,iu), 

and for all u G fl^, there exists a compact absorbing set at time in H for 
ip\t,uj). 

Proof. We have proved in Lemma [241 that the ball B{0, Rslu) is an ab- 
sorbing set at time in V, which is compact in H. Hence, defining K{uj) : = 

|m G y : WuW^ < -R3 + ||2^°(0, ti^) |||,| concludes the proofs. ■ 

Theorem 26 For every value of the parameter A G M, the random dynamical 
system ip^ associated to the equation ^ has a unique global random attractor 
Axiuj). 

Proof. Using Theorem [17] and the existence of a compact absorbing set 
in H, we have the existence of a random attractor Ax{uj) which is forward 
invariant. ■ 
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We can now apply Theorem 2 from [S]. The statement of this theorem 
is composed of two parts, the first one devoted to the convergence of the 
random attractor to the deterministic one as the intensity of the noise goes 
to zero; the second one to the upper semicontinuity of the random attractor 
when the parameter of the noise varies with continuity to some non-zero 
value. We apply the second part. The assumptions of the second part are: i) 
the existence of the random attractor for every fixed value of the parameter, 
ii) the P-a.s. continuous dependence of trajectories on the parameter, in any 
fixed finite interval of time, uniformly in the initial conditions taken from 
any fixed non-random bounded set. Both assumptions have been proved in 
the previous sections. Thus we get the following final result. 

Theorem 27 Let Ax {uj) be the random attractor associated with equation 
iQ, then there is upper semicontinuous convergence of A\{uj) to Ao{t^) as 
X 0: 

lim dn (Ax (u) , An (u)) = with P — a.s. . 
4.3 Random squeezing property 

In this section, we are going to establish that the random attractor of the 
random dynamical system associated with equation has a finite Haus- 
dorff dimension (notice here that for simplicity of notation, we dropped the 
superscript Xin ip^. Let u and v be two solutions of the associated equation 
iQ, then the difference u — v is solution of 

d(u — f ) ^ ^ ^ 

h lyAiu — v) = —B\{u, u — v) — B\{u — v,v). (38) 

(JjL 

Lemma 28 Let U be the orthogonal projection on the first n eigenvectors of 
the operator A. Then, 

^ I R^{e,uj)ds 

\Il{tfi{t,u)uQ- tfi{t,uj)vo)\fj <\uq- volfjc Jo (39) 



I (/ - n) ((^(t, Uj)Uq - ip{t, Uj)Vo) \jj < 



( 



~kn + iut 



[Ri{9s{u))f + ^Ri{9s{uj))ds 



(40) 
\ 
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where C^, is the constant in inequality (jH]), and C is the constant in the 
inequalities in Lemma\^ for all t > and all Uq,Vq G A{u), where Ri = 
Ri ~'^^ given in Lemma l^ say for a given fixed a = Oq := 2a^{u) (see ^E^))■ 

Proof. We multiply equation (138|1 hj u — v, and use inequality ([6]) to 
obtain 

+ ^11(^-^^)111 < I < {Bxiu-v,v)),iu-^ > I 

<^\\{^-v)\\l + ^\\Mu-v,v)\\y, 

T^uf— ~\||2 I C** |~ ~|2 |~|2 

< I^\\{u-V)\\y+—\U-V\f^\v\^. 

Using Gronwall lemma, we obtain 

\u{t)-v{t)\\ < \u{t,)-v{t,)\\e^O^^'^\V' 

Now, using the invariance of the attractor, if we take uq^vq G A{uj) then 
u{t),v{t) G A{6tuj), and therefore by Lemma [21] and Lemma [221 we have 

\ip{t,u;)uo-ipit,iu)vo\l < \uo-vo\%e^^oR^i<^s^)'i^ (41) 

for all f > and all uo,vo G ^(u;). 

Recall that 11 is a projection on the n-dimensional subspace of eigenvec- 
tors of the operator A, we have 

|n {ip{t, Lj)uo - (p{t, Lj)vo) || < \ip{t, Lj)uo - (p{t, Uj)vo\%. 

Let Q := / — n, it commutes with A but not with B\. 

Let us apply the operator Q to the equation ([38l) , then using Lemma [2] and 
the Poincare inequality we get 
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+ 1 < Q{Bxiu-v,v)),Qiu~v) > I 

< (^\Bx(u,u-v)\j^ + \Bx(u-v,v)\jj^ \Q(u-v)\j^ 

< -^=\U-V\h {\\u)\\y + \\v)\\y) \\Q{U-V)\\y 

Now, using the Poincare inequality on the left side of the above inequality 
we get that 



— |g(ix -v)\% + ukn+i\Q{u -v)\l< - v\% {\u\l + \v\l) . 

Hence, 

\Q{u-v){t)\l< \Q{u-v){to)\l-kn+ii^ f \Q{u-v){s)\%ds 

J to 

\{u-v){s)\l{\u{s)\l + \v{s)\l) ds 



C2 



Using Gronwall lemma we get 



to 



\Q{u-v){t)\\ < |g(5r-?;)(to)|2.e-^"+i'^(*-*o) 

+ ^ /*e-'=--(*-^)|(iI-^y)(.)H (1^1(5)11 + \v{s)\l) ds. 



to 



Let us take to = and uq^vq G A{uj) , then we have, thanks to Lemma [2T] 
and Lemma [221 that 



\Q{u-v){t)\l<\u,-Vo\\e-'-^-^'+-^ I e-^-^^'^^'-^^R,{9M)\{u~v){s)\lds. 

kn+lV ' 
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We use ( HTi) in the above inequality to obtain 

\Q{u-v){t)\l < \uo-vo\l (e-''-+'^' + re-^"+i"(*-^)i?i(^,(a;))e^o 

V K+ii^ V /Jo 

On the other hand using Cauchy-Schwarz inequahty and the fact that ^/x < 
for all x > 0, 



1/2 / rt 



[Ri{es{uj))f ds 



1/2 



< 



1 \i/2 / [/2i(^,(a;))]^rfs 



Combining all the above estimates we get that 

/ 



\Q{u-v){t)\l < \u„-v„\l 



-k„+ivt _|_ 



\ 



3/2 



[R,[e,{u;))f + ^R,{e,{u))ds^ 



4.4 Finite dimensionality of the random attractor 

In order to be able to apply Theorem [201 we need to show that EICt-c^uj)) < 
oo, where 

Cniuj):=[R,iuj)f + ^R,{u), (42) 

is the exponent in equation fj40l) of the squeezing Lemma [281 This is because 
C-H^uj) plays, in our case, the role of 6*5(0;) in Theorem [201 Here Ri{uj) = 
given in Lemma [221 (see also Lemma [25]) . say for a given fixed 
a = ao := 2a^,(z/) (see ( 125]) ). In order to get the finite expectation of C-h, i.e. 
E^Cfii^^)) < 00, we need to estimate the moments of the radii Ri. 
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Proposition 29 Let > 6e the constant in the inequality 70 = uko, 
and aQ = 2«^,(z/) (see Then the stationary process % that solves of the 

equation 

dzt = - (a + ao) z^dt + dWt 



satisfies 



E 



for all s < T < 0, where C2 = E 



< 00. 



Proof. By Young's inequality we liave 



C* W^tWy 



\\zt\\y = 2A/7o(r - s) 



T — S] 



< loir - s) + 



^2 ||~||2 



V 



470 (r - 8) 



therefore, 



E 



< e^°^^-'^E 



g (t-s) J a 47Q 



g (t-s) Js 470 



Thus, it is sufficient to show that 

E 

Thanks to Jensen inequality we have 
E 



(r — s) J s 



g (t-s) Js 47Q 













e "'to 


/ s 





Since Zt is a stationary process then it follows that 



E 



e ■'to 



E 



< 00 , 



thanks to (1281) . The proof is complete. 
Let us recall that 



J —00 
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where 



m :=i^i?(t)r+|^|i-(t)r, 



and that Ri{uj) = Rf~'^\uj), which, as it has been remarked in section 
depend on the parameter a. Hereafter we choose a = ckq = 2a=K(z/) (see f l28l) 
and Propositon [29]) . Therefore, the relevant results of section 142) are valid 
for this choice of a. 

Lemma 30 Let 70 and ao be as in Proposition Then 

E{{R^{uo)f) < 00, 
and by the Cauchy-Schwarz inequality 

ERiiuj) < 00. 

Consequently, 

EiCn)<oo, (43) 
where Cji{uj) is given in equation (J^- 



43 



Proof. First we observe that using successively Jensen inequality, Fubini 
Theorem and Holder inequality yield 



1+ f{s)e 



<2 + E 



f{s)e 



2 / ds 



ds 



\ 



ro 2 I { 2CJz(r 

<2 + E f{sfe 

( 



2 



dr 



ds 



<2+ / E 



2 / [2C4z{r)\\-^\dr 



ds 



<2+ (Efisr) 



4x1/2 



^^aJ' ^2C4z{r)\\-^] dr^ 



1/2 



v 



ds. 



The process z is stationary, thus 

E{f{sr) = E{fm. 

All moment of a Gaussian random variable are finite, hence 

E{f{0Y) < oo 

Now, using the preceding estimates and Proposition [29] we conclude that 



44 



E{R,{u))'<2+{Ef{0)') 



4x1/2 



<2 + {Ef{oyy/' f 

J —c 



Ee 



(8CJ|Hr)|| -2kou) dr 



1/2 



ds 



( 



Ee-Js 



^CA\zir)\\dr 



1/2 



ds 



Since 70 = vk^, then 



" {kov - —)s 



2 ' ds. 



E{{R,{uj)Y) < 00, 

and the proof is complete. ■ 

As a consequence, we have the following theorem 

Theorem 31 Let Ki, K2 and be the absolute constant stated in Theorem 
I2U Let n be large enough such that 



n+lj 



\3/2 



< Ki, and kn+iv > K2E (Cn) , 



where C-^ is given in fH2]) . for a = ao = (see / f^) ). Then, P-a.s. the 

random attractor Ax{uj) of the random dynamical system Lp^ associated with 
equation (Q has finite Hausdorff dimension which is less than K^nlnn. 



Proof. The proof follows from applying Theorem [20] for /i 



(l/kn + l) 



3/2 



and 



6 = kn+iJy. Then by virtue of (H3ll all the assumptions of Theorem [20 are 
satisfied. Hence, we get that P-a.s. the random attractor ^a(^) of the ran- 
dom dynamical system cp^ associated with equation (Q has finite Hausdorff 
dimension which is less than K^n In n. This completes the proof. ■ 
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